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I.  INTRODUCTION 


Applied  Research  Laboratories,  The  University  ot  Texas  at  Austin  (ARL:UT), 
has  been  investigating  both  theoretically  and  experimentally  the  basic  properties  of 
thermoacoustic  sources.  The  thermoacoustic  process  employs  direct  heating  of  the 
acoustic  medium  to  produce  a  controlled  local  thermal  expansion  which,  in  turn, 
generates  sound  waves.  The  thermal  energy  is  delivered  by  a  laser  beam  to  the 
water  through  the  free  surface  without  physical  contact.  The  efficiency  of  the 
thermoacoustic  conversion  process  is  an  important  issue.  An  analysis  was  made  of 
the  upperbound  and  it  was  found  to  be  constrained  by  the  physical  properties  of  the 
medium  and  the  energy  output  of  the  thermal  source.1 

The  thermoacoustic  process  is  only  one  of  several  optoacoustic  processes  for 
converting  optical  energy  into  acoustic  energy.  It  is  often  referred  to  as  the  linear 
process.  It  is  linear  in  the  sense  that,  over  a  large  range  of  temperatures,  the  change 
in  density  of  the  water  is  linearly  proportional  to  the  optical  energy  input,  at  least  for 
small  input  energy  densities.  The  general  concept  is  illustrated  in  Fig.  1.  Other 
processes  by  which  optical  energy  is  converted  to  sound  are  grouped  together  under 
the  title  of  nonlinear  processes.  The  current  knowledge  regarding  nonlinear 
processes  was  thoroughly  reviewed  by  Lyamshev  and  Naugol'nikh.2  Nonlinear 
processes,  in  general,  have  been  found  experimental ly  to  be  more  efficient  but  they 
are  not  well  understood. 

The  anatomy  of  the  optoacoustic  processes  is  illustrated  in  Fig.  2.  The  main 
processes  include:  weak  and  strong  evaporation,  optical  breakdown  of  the  vapor, 
and  optical  breakdown  of  the  liquid.  They  are  mainly  governed  by  the  intensity  and 
energy  density  of  the  laser  pulse.  Evaporation  processes  are  mainly  governed  by 
energy  density.  Weak  evaporation  refers  to  the  evaporation  of  a  thin  surface  layer  of 
the  liquid,  while  strong  evaporation  refers  to  explosive  bubbling  from  the  surface  and 
from  the  interior  of  the  liquid. 

Evaporation  processes  can  take  place  only  when  the  energy  absorbed  is  equal 
to  or  greater  than  the  specific  internal  energy  required  for  the  phase  change.  The 
absorbed  energy  is  directly  governed  by  the  delivered  surface  energy  density  and  by 
the  absorption  coefficient.  For  1.06  pm  laser  radiation,  the  absorption  coefficient  is 
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normally  about3  15  m"1,  which  gives  a  surface  energy  density  threshold  of 
approximately  2  x  107  J/m2. 

Above  certain  intensity  and  energy  density  thresholds,  optical  breakdown 
causes  the  physical  properties  of  the  medium  to  change.  The  intensity  threshold  for 
breakdown  of  water  is  strongly  dependent  on  impurities  such  as  suspended 
particles.  The  breakdown  causes  the  water  to  become  opaque.  This  is  equivalent  to 
a  sudden  increase  in  the  absorption  coefficient^  and  a  consequent  reduction  in  the 
evaporation  threshold. 

For  10.6  pm  laser  radiation  on  water,  optical  breakdown  of  the  ejected  vapor 
has  been  observed3  above  a  threshold  of  2  x  1012  W/m2.  Of  course,  there  must 
also  be  sufficient  energy  density  to  produce  the  vapor  in  the  first  place.  Therefore 
there  is  an  energy  density  threshold  as  well. 

Sound  generation  is  by  a  combination  of  two  main  physical  processes,  the 
expansion  of  the  material  within  the  liquid  due  to  physical  or  chemical  changes  and 
the  recoil  pressure  at  the  surface  from  the  ejected  vapor.  There  have  been  a  number 
of  publications  on  this  topic.  A  large  proportion  of  the  reported  works  have  been  for 
cases  where  most  of  the  energy  is  spread  over  several  megahertz.6-7  The  results 
reported  by  Hickman8  were  in  the  range  of  sonar  frequencies  but  the  experimental 
conditions  were  imprecise.  The  most  extensive  set  of  results  were  reported  by 
Maccabee  and  Bell.5*9*  10,11  Their  data  covered  a  wide  range  of  operating 
conditions.  Their  results  showed  significant  finite  amplitude  distortion  and  extra 
attenuation  of  the  acoustic  pressure.  Maccabee  and  Bell  used  a  "rocket"  model  to 
predict  the  blast  pressure  and  linear  superposition  to  construct  the  acoustic  pressure. 
They  found  that  the  model  predictions  were  not  consistent  with  the  experimental 
results.1 1  One  likely  contributor  to  the  discrepancy  is  the  nonlinear  propagation  of 
the  acoustic  output,  because  they  used  laser  pulses  of  very  high  energy  densities 
and  intensities. 


II.  THEORY 


The  optoacoustic  process  being  considered  here  is  specifically  th..:  cf  •  •  • 
breakdown,  followed  by  a  phase  change  from  liquid  to  gas  and  jr o  w; 'e  .  .< 
expansion  of  the  gas.  The  initial  dielectric  breakdown  and  phase  c  ' ■  c * 

expected  to  significantly  contribute  to  the  sound  generation  process  but  v  oy  w 
consume  a  certain  amount  of  energy.  This  may  be  considered  as  an  overhead 
requirement.  The  explosive  expansion  is  believed  to  be  the  mam  mechan.sm  cf  sound 
generation  It  will  be  modeled  as  a  blast  reaction.  The  blast  model  used  is  based  on 
the  work  of  A.  N.  Pirri1^'1^  who  successfully  modeled  the  momentum  transfer  from  a 
high  powered  laser  blast  to  the  surface  of  a  solid.  The  model  is  general  enough  to  be 
adapted  to  give  a  first  approximation  of  the  acoustic  signals  produced  by  a  laser 
induced  blast  at  the  water  surface.  Our  model  is  based  on  the  more  complicated  of  the 
two  Pirri  models.  It  is  a  two-stage  model.  It  has  a  one  dimensional  stage  where  the 
blast  front  is  assumed  to  rise  vertically  from  the  surface.  This  stage  is  essentially  the 
"rocket  model"  used  by  Maccabee  and  Bell.1 1  The  second  stage  is  a  two-dimensional 
spreading  regime  where  the  blast  expands  outward  over  the  surface.  The  model  is 
illustrated  in  Fig.  3.  The  sound  is  assumed  to  propagate  linearly  and  the  underwater 
acoustic  pressure  signal  is  constructed  from  elementary  constituents  by  linear 
superposition.  In  a  theoretical  study  without  experimental  verification,  Wu14  followed 
this  approach  and  obtained  an  estimate  for  the  acoustic  pressure.  The  deformation  of 
the  water  surface  is  assumed  to  be  negligible.  The  linear  propagation  and  negligible 
deformation  assumptions  are  some  of  the  main  shortcomings  of  th.s  model 
Nevertheless,  under  operating  conditions  where  these  assumptions  are  approximately 
valid,  the  model  is  expected  to  produce  useful  results.  With  further  theoretical  analysis 
supported  by  experimental  measurements,  more  realistic  models  will  evolve 

A.  Basic  Blast  Theory 

The  nonlinear  optoacoustic  reaction  is  modeled  as  an  explosive  blast  generate,; 
by  the  laser  radiation  at  the  water  surface.  The  blast  is  assumed  to  be  gas  driven  For 
simplicity,  let  us  assume  that  the  gas  within  the  blast  is  homogeneous  Given  that  the 
rate  of  energy  input  is  sufficient  to  cause  the  blast  wave  to  be  supersonic,  the 
relationship  between  the  internal  energy  Eg,  the  pressure  Pg,  and  volume  Vg  of  the 
blast  gas  is  given  by  the  adiabatic  gas  equation, 
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where  y  is  the  ratio  of  specific  heats  of  the  blast  gas. 


The  movement  of  the  blast  front  is  governed  by  conservation  of  momentum 
considerations.  The  difference  Pqw  between  the  blast  pressure  Pg  and  the  ambient 
pressure  P0  must  be  equal  to  the  rate  of  change  of  momentum  per  unit  area  at  the 
blast  front;  therefore, 

PDW  =  pg  -  po  =  (du/dt)2  po/(y  +  1 )  ,  (2) 

where  p0  is  the  density  of  the  surrounding  fluid  and  u  is  the  distance  along  the  outward 
normal  at  the  boundary  of  the  blast. 

Assuming  isentropic  flow,  the  coupling  coefficient  Cc  between  the  blast  pressure 
Pg  and  the  pressure  Ps  at  the  water  surface,  quoting  from  Ref.  13,  is  given  by  the 
relationship, 

Cc-Ps'Pow  =  !<r*i>'2rf2i'/<lM)  •  (3) 

These  three  relationships  form  the  basis  of  the  blast  theory  used  by  Pirri  to 
estimate  momentum  transfer  from  a  high  powered  laser  pulse  to  a  solid  surface.13 
Pirri’s  model  has  been  modified  to  predict  the  acoustic  waveform  produced  by  a  high 
powered  laser  pulse  incident  on  a  water  surface. 

B.  Construction  of  the  Acoustic  Signal 

The  acoustic  signal  produced  by  a  high  powered  laser  beam  striking  the  surface 
of  the  water  will  be  constructed  by  superposition.  The  laser  beam  is  assumed  to  have 
a  uniform,  circular  cross  section.  For  simplicity,  the  acoustic  signal  is  observed  at  a 
distance  that  is  sufficiently  large  to  be  outside  the  nearfield  of  the  source  to  allow  the 
use  of  farfield  approximations. 


The  blast  model  program  was  implemented  as  a  modular  FORTRAN  program  that 
analyzes  the  received  signal  in  three  stages.  First,  the  blast  pressure  and  area  of  the 


blast  produced  by  the  laser  are  computed  as  functions  of  time.  Second,  the  program 
computes  the  time  derivative  of  the  downward  pressure  at  the  surface  due  to  the  blast 
using  the  blast  pressure  and  area  histories.  Third,  the  acoustic  signal  at  a  given 
distance  and  angle  is  computed  from  the  pressure  derivative  using  Green's  function 
and  linear  superposition.  Spreading  and  absorption  losses  are  taken  into  account. 
Movement  of  the  water  surface  is  neglected.  Pirri  initially  used  a  spherical  blast 
model,1  2  but  he  later  replaced  it  with  a  two-stage  model,13  shown  in  Fig.  3.  The 
current  model  developed  in  this  report  is  based  on  the  latter.  The  laser  pulse  is 
assumed  to  be  of  constant  power  and  beamwidth  with  a  finite  duration  Tp  The  blast 
model  was  divided  into  two  stages  as  shown  in  Fig.  3.  In  the  first  stage,  the  blast  is  a 
thin  layer  of  vapor  spreading  upwards  from  the  surface.  Since  initially  the  thickness  is 
small  compared  to  the  width  of  the  blast  area,  the  expansion  is  mostly  upwards  and 
one  dimensional,  with  negligible  edge  effects.  In  the  one-dimensional  stage  the 
velocity  of  the  blast  front  velocity  du/dt  is  denoted  by  Vq^.  Applying  Eqs.  (1)  and  (2) 
for  the  one-dimensional  case  and  assuming  that  the  laser  intensity  and  beamwidth  are 
constants,  it  can  be  shown  that  Vpyy  is  also  a  constant  The  expression  for  Vq^, 
originally  derived  by  Raizer,15  is 

vDw  =  (2<i2-')|o/Po]v3  <«) 

where  l0  is  the  laser  intensity  in  W/m2.  Substituting  for  du/dt  in  Eq.  (2)  by  the 
expression  for  Vq^/  from  Eq.  (4),  the  excess  pressure  within  the  blast  volume  is  given 
by 


PDW  =  PovDW2  ;  (Y+ 1)  (5) 

Applying  the  coupling  relationship  given  in  Eq  (3),  the  pressure  ps1Q  at  the  water 
surface  in  the  one-dimensional  stage  is  given  by 

PslD  =  PDWcc 

If  the  laser  power  is  turned  off  while  the  blast  is  still  one-dimensionat.  the  pressure  will 
decay  as  a  planar  blast1®  where  the  pressure  ps(t)  is  given  by 

PsM  =  PslD  ( V*p  y2'3  •  for(Tp<t<  t2D>  <7> 


- -  -  ■ 
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where  tp  is  the  duration  of  the  laser  pulse  and  x2d  is  the  duration  of  the 
one-dimensional  stage.  When  the  blast  layer  thickness  approaches  the  width  of  the 
blast  area,  edge  effects  are  no  longer  negligible  and  the  pressure  to  expand  outwards 
can  no  longer  be  ignored.  At  this  point,  the  blast  enters  its  second  stage.  The  second 
stage  is  a  two-dimensional  spreading  regime  where  the  blast  expands  outward  over 
the  surface.  The  blast  geometry  is  assumed  to  change  from  one-dimensional  to 
two-dimensional  when  the  height  of  the  one-dimensional  blast  equals  its  diameter. 
This  transition  occurs  at  time  t20,  given  by16 

A 

x2D=  (1/V)  dr  ,  (8) 

o 

where  Ds  is  the  diameter  of  the  one-dimensional  blast  area,  also  equal  to  the  diameter 
of  the  laser  illuminated  area  on  the  water  surface,  r  is  the  displacement  normal  to  the 
wavefront,  and  V  is  the  velocity  of  the  blast  wave  identically  equal  to  du/dt  in  Eq.  (2).  If 
tp  >  t2d.  then  the  velocity  V  of  the  blast  front  is  equal  to  the  constant  and  x2d  is 
given  by 

t2d  =  ds/vDW  ,  (9) 

If  tp  <  x2d  then  the  velocity  V  of  the  blast  front,  after  the  cessation  of  the  laser  power 
input,  is  given  by  the  planar  decaying  blast  theory.16  Quoting  from  Ref.  13,  t2d  is 
then  given  by 

'2D  -  l2/(3  Vow312  'p1/2)J  [Ds3/2  -  (VDW  tp)3/2J  +  Tp  ;  for  (Tp  £  t  <  t2D)  .  (10) 

After  time  x2d  the  blast  is  assumed16  to  decay  according  to  cylindrical  blast  theory 
irrespective  of  any  additional  laser  input  energy.  The  equation  for  the  pressure  is 
given  by 

PS  *  Psl  D  (x2p 1  t2D)2/3  (x2D  7 1)  *  (11) 

where  t2p=tp  if  tp<t2o  or  x2p=x2D-  ^Xp-T2D-  An  illustration  of  the  surface 
pressure  as  a  function  of  time  is  given  in  Fig.  4, 


Using  the  expression  for  a  pressure  source,  the  acoustic  pressure  p|_  at  a  point 
(x0,  y0,  z0)  in  Cartesian  coordinates  in  a  lossless  medium  is  given  by 

PL  *  J'J(dps(t-R'c)/9t)  /  (2rccR)  dy  dx  .  (12) 

It  is  assumed  that  the  blast  is  cylindrical^  symmetric  with  a  radius  r  which  varies  with 
time  t.  R  is  the  distance  from  the  point  (x,  y,  0)  on  the  surface  to  the  hydrophone  at  (x0, 
y0,  z0).  For  simplicity,  and  without  loss  of  generality,  let  the  blast  area  be  centered 
about  the  origin  and  y0  =  0.  Therefore,  R  is  given  by 

R  =  ^((x-x0)2  +  y2  +  (z-z0)2)  •  (13) 


Then,  the  distance  R0  from  the  center  of  the  blast  area  to  the  hydrophone  is  simply 


When  the  hydrophone  is  in  the  farfield,  R  may  be  approximated  by 


R  -  V«x  '  xo)2  ♦  ( z  *  zo)2) 


(14) 


which  is  independent  of  y.  Assuming  that  the  blast  gases  are  well  mixed,  the  surface 
pressure  and  its  rate  of  change  will  be  uniform  within  the  blast  area.  Assuming  also 
that  the  pressure  at  the  boundary  of  the  blast  area  decays  to  zero  over  a  finite  transition 
region  or  skirt  of  width  h,  the  expansion  of  the  integral  with  respect  to  y  for  any  time  t 
may  be  expressed  as 


(dps/t)t)  dy 


(y1-(hr/2y1)) 

(dps/dt)dy  +  2(hr/y1)(ps/r)Vs 


~(yi  -(hr/2y  i )) 


(y-|  +(hr/2y1 )) 

+  2  J  (dpg/dt)  ((y-CV!  -(h r/2yi  ))/(hr/yi ))  dy  ,  (15) 

(yi(hr/2y1)) 

where  Vs  is  the  surface  velocity  of  the  blast  boundary,  and  the  limit  y-|  is  a  function  of 
the  value  of  x  and  the  radius  r  of  the  blast  area, 


yi  =  V(f2  -  x2>  • 


Figure  5  shows  a  model  of  the  surface  pressure  distribution.  The  first  term  in  Eq.  (16) 
is  due  to  the  rate  of  change  of  blast  pressure  within  the  homogeneous  blast  area,  the 
second  term  is  the  contribution  from  the  rate  of  change  of  pressure  at  a  point  within  the 
skirt  due  to  the  movement  of  the  blast  boundary,  and  the  last  term  is  due  to  the 
pressure  change  in  the  skirt  of  the  blast.  Equation  (15)  simplifies  to 


(aps/dt)  dy  *  (3ps/3t)  y1  +  rp^y,  , 

J 


which  is  independent  of  h,  the  width  of  the  skirt  of  the  blast  area.  The  acoustic  signal  at 
a  given  distance  and  angle  is  computed  from  the  downward  surface  pressure  by  taking 
into  account  the  spreading.  Since  the  time  delay  term  R/c  is  a  function  of  x,  the  radius  r 
of  the  blast  is  also  dependent  on  x  and  the  relationship  between  (3ps(t-R/c)/3t)  and  x 
becomes  rather  complicated  particularly  when  the  blast  velocity  is  supersonic; 
therefore,  the  remaining  integral  with  respect  to  x  is  carried  out  numerically.  The  blast 
surface  pressure,  its  time  derivative,  and  the  blast  radius  as  functions  of  time  are 
generated  as  spline  curves  to  ensure  continuity  of  the  first  derivative. 


After  the  acoustic  signal  for  a  lossless  medium  has  been  constructed,  the 
absorption  losses  are  calculated  using  the  frequency  dependent  absorption  coefficient 
a  in  dB/m  from  Schulkin  and  Marsh,17 

a  =  (  ASfjf2  )/( fT2  +  f2  )  +  Bf2/fT  ,  (18) 


FIGURE  5 

MODEL  OF  THE  SURFACE  PRESSURE  DISTRIBUTION 


f  =  frequency  in  kHz, 

fT  -  21.9  x  io(6-1-52/(T+273)).  (19) 

A  =  0.00002034, 

B  -  0.00002931, 

S  ■  salinity  in  parts  per  thousand,  and 
T  *  temperature  in  degrees  Celsius. 

The  absorption  coefficient  is  then  applied  to  the  signal  in  the  frequency  domain  aftc 
taking  its  Fourier  transform.  The  final  signal  is  then  reconstructed  from  its  transform. 

C.  Theoretical  Results 

An  example  of  the  theoretical  estimate  of  the  acoustic  signal  waveform  and  its 
energy  density  spectrum  is  shown  in  Figs.  6(a)  and  (b),  respectively.  Referring  to  the 
geometry  in  Fig.  5,  the  range  from  the  source  to  the  point  of  observation  R0  is  1 .2  m 
and  the  ray  path  makes  an  angle  with  the  downward  vertical  0O  of  45°.  Since  the  blast 
very  rapidly  enters  the  second  stage  and  expands  outward  supersonically,  the 
contribution  from  the  edge  of  the  expanding  and  decaying  second  stage  of  the  blast 
arrives  before  the  main  contribution  from  the  more  impulsive  first  stage,  thus  giving  the 
acoustic  pulse  a  gradual  buildup  in  pressure  before  the  arrival  of  the  main  pressure 
pulse.  The  main  pressure  pulse  is  very  short,  and  is  followed  by  a  rarefaction  phase 
to  complete  the  cycle.  The  energy  density  spectrum  plot  shows  the  energy  mainly 
concentrated  below  200  kHz.  Below  50  kHz  the  spectrum  falls  off  due  to  the  transfer 
function  of  the  optoacoustic  process,  which  decreases  with  decreasing  frequency. 
Above  100  kHz  the  spectrum  abruptly  falls  off  due  to  diffraction  and,  to  a  lesser  extent, 
absorption  losses.  The  absorption  loss  increases  with  frequency  and  range.  The 
diffraction  effect  becomes  significant  when  the  resolved  physical  size  of  the 
optoacoustic  source  in  the  direction  of  the  ray  path  from  the  source  to  the  point  of 
observation  exceeds  the  acoustic  wavelength.  The  model  assumes  that  the 
optoacoustic  source  is  infinitely  thin.  Therefore,  any  diffraction  effects  must  come  from 
the  re  solved  dimensions  of  the  blast  area  in  the  direction  of  the  ray  path;  the  resolved 
length  of  the  source  in  this  direction  is  approximately  the  diameter  of  the  blast 
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FIGURE  6  (a) 

AN  EXAMPLE  OF  THE  THEORETICALLV  PREDICTED 
WAVEFORM  FOR  THE  FOLLOWING  CONDITIONS: 

R0=  1.2m 
80=  45- 

laser  pulse  energy  =  0.1 3  J 

laser  beam  area  =  8  mm^ 

laser  pulse  duration  =  3.5  jis 

ratio  of  specific  heats  of  blast  gas  =1.2 
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FIGURE  6(b) 

THE  ACOUSTIC  ENERGY  DENSITY  SPECTRUM  OF  THE 
PRESSURE  WAYEFORM  IN  FIG.  6(a) 
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area  multiplied  by  sin  (90),  where  0O  is  the  angle  between  the  downward  vertical 
and  the  hydrophone.  This  explanation  was  found  to  be  consistent  with  the  observed 
upper  cutoff  frequency. 

D.  Comparison  of  Theoretical  Predictions  with  the  Results  of  Maccabee  and  Bell 

Maccabee  and  Bell ^  ^  presented  a  number  of  experimental  results  of  peak 
acoustic  pressure  as  a  function  of  laser  pulse  energy  at  the  109th  Meeting  of  the 
Acoustical  Society  of  America.  One  set  of  their  data  satisfies  the  farfield  criterion  for 
which  the  present  model  is  applicable.  The  comparison  between  the  peak  pressures 
predicted  by  our  model  and  their  experimental  data  is  shown  in  Fig.  7.  It  can  be  seen 
that  there  is  some  agreement  between  them,  although  the  theoretical  and  experimental 
waveforms,  which  are  not  shown  here,  were  quite  different.  The  difference  in  the 
waveforms  was  expected  because  our  model  assumes  only  linear  propagation  while 
the  experimental  data  appeared  to  be  shock  limited.  Nevertheless,  within  the  band  of 
the  signal  frequencies  and  over  such  a  short  distance,  the  extra  attenuation  due  to 
nonlinear  propagation  was  not  expected  to  be  very  significant  and  therefore  the 
comparison  of  the  peak  pressure  of  the  leading  edge  of  the  pulse  is  expected  to  be 
valid.  The  results  show  that  at  the  lower  energy  levels,  the  theory  appears  to  over¬ 
estimate  the  peak  pressures.  The  way  the  theory  and  experimental  results  diverge 
suggest  that  a  certain  amount  of  the  laser  pulse  energy  may  not  be  directly  contributing 
to  the  blast.  This  is  not  inconsistent  with  the  current  understanding  of  nonlinear 
optoacoustic  processes  in  which  the  medium  has  to  undergo  a  local  change  of  state, 
such  as  a  dielectric  breakdown,  before  the  blast  can  take  place.  The  change  of  state 
will  consume  a  finite  amount  of  the  laser  energy  leaving  the  remainder  to  fuel  the 
optoacoustic  generation  process  which  is  modeled  as  a  blast.  Therefore,  a  constant 
laser  energy  density  overhead  should  be  expected.  From  Fig.  7,  it  is  estimated  that 
the  energy  density  overhead  is  about  3.3  x  102  J/m2. 
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FIGURE  7 

COMPARISON  OF  THEORETICAL  MODEL  PREDICTIONS 
OF  PEAK  PRESSURE  WITH  THE  EXPERIMENTAL  DATA 
FROM  MACCABEE  AND  BELL 

THE  EXPERIMENTAL  CONDITIONS  WERE: 

Rq=  0.5  m 

e0=  °° 

LASERBEAM  AREA  =  0.2  cm2 
LASER  WAVELENGTH  =  1.06  nm 
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III.  EXPERIMENT 


Experiments  were  carried  out  to  verify  the  blast  model  and  to  obtain  a  better 
understanding  of  the  process.  Due  to  the  small  amount  of  time  and  resources 
available  at  this  stage  of  the  contract,  these  investigations  should  be  considered  as 
preliminary.  In  these  experiments,  the  specific  goals  are  (1)  to  measure  the  energy 
density  overhead,  (2)  to  compare  the  predictions  of  the  blast  model  with  the  practical 
results,  and  (3)  to  identify  weaknesses  in  the  blast  model  which  need  to  be  improved. 
From  these  results,  recommendations  for  future  work  will  be  made. 

A.  Apparatus 

The  experimental  arrangement  is  shown  in  Fig.  8.  The  laser  used  is  an  Apollo 
Nd:glass  laser  Model  22/A  with  a  peak  energy  of  22  J  and  a  wavelength  of  1 .06  pm.  It 
is  capable  of  producing  an  unmodulated  pulse  of  up  to  1  ms  and  a  spot  size  of 
approximately  8  mm2.  A  short  laser  pulse  is  achieved  by  putting  a  Pockels  cell  and 
polarizer  in  the  laser  cavity  and  using  a  Q-switch  driver  to  modulate  the  Pockels  cell. 
Each  pulse  lasted  less  than  3.5  ps  and  a  train  of  up  to  five  pulses  at  a  maximum 
repetition  rate  of  10,000  pulses  per  second  could  be  delivered.  The  control  circuit 
determined  the  number  of  pulses  and  the  timing  between  pulses.  A  fast  photodetector 
connected  to  a  Nicolet  digital  oscilloscope  monitored  the  laser  pulse.  The  laser  beam 
was  directed  vertically  down  to  the  water  surface  by  a  prism  and  focused  on  the 
surface  by  a  pair  of  lenses.  The  acoustic  signal  was  detected  by  a  hydrophone.  The 
hydrophones  used  were  the  E8-59,  E8-62,  and  H23.  Detailed  descriptions  of  the 
hydrophones  are  given  in  Appendix  A.  The  outputs  of  the  hydrophones  were  amplified 
and  filtered  through  a  30  dB  balanced  amplifier  or  an  HP465A  20  dB  amplifier  and 
followed  by  a  KH3100  bandpass  filter.  The  data  were  recorded  on  a  Nicolet  4096 
digital  oscilloscope. 

B.  Noise  and  Multipath  Suppression 

In  our  experiments,  we  encountered  numerous  noise  interference  problems.  The 
most  persistent  of  these  was  the  electrical  noise  pulse  from  the  Q-switch  driver.  This 
noise  pulse  affected  the  performance  of  the  control  circuit  and  appeared  in  the  output 
signals  of  the  photodetector  and  the  amplifiers.  The  Q-switch  driver  sent  a  sharp 
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negative  voltage  spike  through  the  power  supply  of  the  control  circuit  which 
occasionally  caused  the  circuit  to  go  into  the  hardware  equivalent  of  an  infinite  loop. 
This  problem  was  solved  by  putting  a  power  diode  in  series  with  the  power  supply. 
We  eliminated  the  noise  in  the  photo-detector  by  building  two  spatially  symmetric, 
voltage  balanced  photodetector  circuits.  One  was  open  to  the  laser  light  and  the  other 
was  covered.  The  difference  of  the  two  signals  was  the  detected  optical  pulse  of  the 
laser.  A  description  of  the  circuit  is  given  in  Appendix  B. 

The  acoustic  signal  was  contaminated  by  both  electrical  and  acoustic  noise.  Due 
to  the  time  delay  between  the  firing  of  the  laser  and  the  arrival  of  the  acoustic  signal, 
the  direct  electrical  interference  from  the  laser  system  could  be  gated  out.  However, 
residual  ringing  from  the  laser  interference  persisted  long  enough  to  contaminate  the 
acoustic  signal  as  a  low  frequency  background.  This  was  removed  by  digital  filtering. 

Acoustic  and  electrical  noise  from  all  other  sources  were  random  and  they  were 
reduced  by  taking  ensemble  averages.  An  ensemble  average  of  ten  was  usually 
enough  to  yield  a  clean  signal,  with  a  signal-to-noise  ratio  (S/N)  of  about  40  dB. 

Multipath  interference  problems  were  avoided  by  careful  design  of  the  layout. 
The  first  multipath  to  reach  the  hydrophone  is  the  reverberation  from  the  bar  on  which 
the  hydrophone  is  mounted.  By  putting  the  hydrophone  on  a  small  extension,  the 
multipath  was  sufficiently  delayed  so  it  could  be  gated  out.  The  other  source  of 
reverberation  was  the  walls  of  the  tank.  Fortunately  our  tank  is  large  enough  that  these 
would  arrive  well  after  the  signal. 

c.  Erogsduca 

Before  beginning  the  experiment,  the  laser  power  supply  and  Q-switch  driver 
were  switched  on  and  allowed  approximately  5  minutes  to  stabilize.  Then,  the  pulse 
repetition  sequence  on  the  pulse  controller  was  set  to  the  required  values.  The 
hydrophone  was  put  into  position.  While  running  the  experiment,  the  laser  coolant 
temperature,  the  photodetector’s  output,  and  the  received  acoustic  signal  were 
monitored.  The  E8  hydrophones  were  calibrated  from  100  kHz  to  1  MHz,  and  in  this 
band  their  response  is  known  to  be  flat.  Their  outputs  were  bandpass  filtered  from 
75  kHz  to  500  kHz  to  avoid  aliasing  because  the  data  recording  system  could 


only  sample  at  2  MHz.  The  H23  hydrophone  was  calibrated  from  1  kHz  to  150  kHz 
and  within  this  band  its  response  was  approximately  flat;  its  output  was  bandpass 
filtered  over  the  same  band. 


D.  Experimental  Results 

The  acoustic  signal  energy  density  was  measured  at  a  hydrophone  located 
approximately  1  2  m  directly  below  the  source  (R0  =  1.2  m  and  0O  =  0  in  Fig.  8).  The 
measured  acoustic  energy  density  is  plotted  as  a  function  of  corresponding  laser  pulse 
surface  energy  density  in  Fig.  9.  The  data  were  from  a  set  of  multiple  pulse  firings  of 
the  laser.  In  each  firing,  a  sequence  of  five  pulses  was  delivered  to  the  water  surface. 
The  pulse  repetition  period  was  0.125  ms.  At  laser  energy  levels  below  4  J/cm^,  it 
appears  that  the  data  lie  approximately  on  a  straight  line  irrespective  of  the  position  of 
the  pulse  in  the  sequence,  except  for  the  lowest  data  point  where,  referring  to  Fig.  2, 
the  laser  energy  density  may  be  too  low  for  a  blast  to  fully  develop.  At  higher  energy 
density  levels,  it  appears  that  there  are  some  anomalies  and  that  the  position  of  the 
pulse  in  the  sequence  has  some  influence  on  the  result.  The  second  and  subsequent 
pulses  appear  to  have  a  slightly  higher  energy  conversion  efficiency  than  the  first 
pulse.  It  is  speculated  that  the  after-effects  of  the  first  pulse  caused  the  energy  density 
overhead  requirement  for  the  subsequent  pulses  to  be  reduced.  Fitting  a  straight  line 
through  the  data  from  the  first  pulses  of  each  sequence,  the  optical  energy  density 
overhead  requirement  for  a  single  laser  pulse  was  found  to  be  16  x  10^  J/m^.  This 
result  does  not  agree  with  the  corresponding  result  from  the  data  of  Maccabee  and 
Bell.  A  possible  cause  of  the  discrepancy  is  the  difference  in  the  laser  wavelengths. 
Maccabee  and  Bell  used  a  laser  of  wavelength  10.6  pm  while  we  used  1.06  pm.  At 
the  latter  wavelength,  the  laser  absorption  coefficient  in  water  is  significantly  lower,  and 
therefore  it  is  logical  to  expect  that  more  energy  will  be  dissipated  into  the  interior  of  the 
water  before  the  start  of  a  nonlinear  reaction  at  the  surface. 

The  peak  acoustic  pressure  was  plotted  against  the  peak  laser  power  to 
determine  the  empirical  relationship  between  them.  The  results  are  shown  in  Fig  10 
Unlike  the  energy  density  comparison,  it  appears  that,  at  laser  peak  intensities  greater 
than  lOOkW/m^,  the  second  and  subsequent  pulses  have  lower  peak  acoustic 
pressures  for  the  same  peak  laser  intensity  The  data  point  at  the  lowest  peak  acoustic 
pressure  and  laser  intensity  appears  to  be  anomalous  As  with  the  lowest  data  point  m 
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FIGURE  10 

EXPERIMENTAL  MEASUREMENTS  OF  PEAK  ACOUSTIC 
PRESSURE  AS  A  FUNCTION  OF  PEAK  OPTICAL  INTENSITY, 
FROM  THREE  PULSE  TRAINS  OF  FIVE  PULSES  EACH.  A 
straight  line  fit  is  shown  through  the  first  pulses.  The 
experimental  conditions  were: 

Rg  =  1 .2  m 

e0  =  o- 

laser  beam  area  =  8  mm  2 

laser  pulse  duration  =  3.5  ys 

receiver  bandwidth  extended  from  75  kHz  to  500  kHz 
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Fig.  9,  it  is  speculated  that,  referring  to  Fig.  2,  the  laser  intensity  may  be  too  low  for  a 
blast  to  be  properly  triggered.  Fitting  a  straight  line  through  the  first  pulse  data  points, 
the  results  indicate  that  a  peak  laser  intensity  of  approximately  3  x  109W/m2  is  the 
minimum  peak  intensity  for  a  blast  to  be  properly  triggered. 


In  Fig.  11,  the  acoustic  pulses  from  single  laser  pulses  of  different  energy  levels 
are  compared.  The  results  in  Fig.  1 1  were  obtained  using  the  E8-59  hydrophone, 
covering  the  frequency  band  from  100  kHz  to  500  kHz.  In  Figs.  12(a)  and  (b),  the 
energy  density  spectra  of  the  pulses  are  compared.  The  spectra  in  Fig.  12(a)  were 
generated  directly  from  the  pulses  shown  in  Fig.  11.  Those  of  Fig.  12(b)  were  from 
pulses  of  similar  energy  detected  by  the  H23  hydrophone  covering  the  frequency 
range  from  1  kHz  to  150  kHz.  The  spectra  at  the  two  energy  levels  were  found  to  be 
of  the  same  shape.  A  doubling  of  the  laser  pulse  energy  produced  an  increase  of 
approximately  6  dB  in  the  acoustic  energy  density  spectrum. 
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The  acoustic  signals  at  angle  60  =  0  and  45°  were  also  measured.  Examples  of 
their  energy  density  spectra  are  compared  in  Fig.  13.  The  spectrum  from  the  0O  =  45° 
case  appears  to  cut  off  at  approximately  1 75  kHz,  while  that  of  the  0O  =  0°  case  cuts  off 
at  approximately  450  kHz. 


The  cutoff  in  the  45°  case  is  due  to  the  frequency  dependence  of  the  source 
directivity  function.  At  wavelengths  approaching  the  resolved  width  of  the  blast  in  the 
hydrophone  direction,  the  contributions  from  different  parts  of  the  blast  area  will 
combine  destructively.  The  width  of  the  blast  can  be  estimated  from  the  -3  dB  cutoff 
frequency  and  the  pressure  distribution  across  the  blast  aperture.  The  pressure 
distribution  across  the  blast  area  could  probably  have  been  deduced  from  a  more 
detailed  set  of  acoustic  measurements  but  at  this  stage,  for  simplicity,  a  simple 
distribution  function  was  assumed.  The  derivation  of  the  equations  is  presented  in 
Appendix  C.  Assuming  an  even  distribution  of  the  surface  pressure  gives  a  width  of 
1.08  mm,  which  is  an  area  of  0.92  mm2.  Assuming  a  Gaussian  distribution  gives  a 
width  of  0.72  mm,  which  is  an  area  of  0.41  mm2.  The  area  we  determined  using  laser 
burns  on  photographic  paper  is  8  mm2.  If  this  model  is  valid,  then  the  acoustic  results 
would  suggest  that  the  laser  energy  is  mainly  concentrated  in  a  much  smaller  area 
than  the  burn  mark  would  indicate.  Although  the  burn  method  of  determining  the  blast 
area  is  somewhat  inaccurate  since  it  does  not  linearly  represent  the  spatial  distribution 
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FIGURE  1  1 

AN  EXAMPLE  OF  THE  EXPERIMENTALLY  MEASURED  ACOUSTIC 
WAVEFORMS  USING  THE  E8-59  HYDROPHONE  AT: 

Rg=  1.2m 

0O-O° 

laserpuise  energy  =  0.22  J  and  0  41  J 

laser  supply  Yoltage  -  4.9  k  Yand  5  9  kY,  respectively 

laser  beam  area  *  8  mm2 


laser  pulse  duration  =  3.5 
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FIGURE  12  (a) 

THE  ENERGY  DENSITY  SPECTRA  OF  THE  ACOUSTIC  PULSES 

IN  FIG  1  1 
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FIGURE  12  (b) 

THE  ENERGY  DENSITY  SPECTRA  OF  ACOUSTIC  PULSES 
DETECTED  BY  THE  H23  HYDROPHONE  OF  THE  SAME 
LASER  ENERGY  LEVELS  AS  IN  FIG.  11 
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FIGURE  13 


COMPARISON  OF  THE  SPECTRA  FROM  THE  0O  =  45 
AND  THE  0O  =0°  CASES 
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of  the  laser  energy,  it  is  difficult  to  believe  that  it  could  be  in  error  by  an  order  of 
magnitude.  A  possible  cause  for  the  discrepancy  may  be  the  deformation  of  the  water 
surface,  which  is  ignored  in  the  current  model.  In  the  published  literature,  small, 
approximately  spherical  depressions  in  the  water  surface  have  been  reported.®  The 
optoacoustic  source  distribution  is  expected  to  conform  to  the  depression,  which  will 
have  the  effect  of  raising  the  diffraction  limited  upper  cutoff  frequency  in  the  90  *  0 
directions.  The  cutoff  at  0°  may  be  due  to  either  the  thickness  of  the  layer  of  effective 
sound  sources  or  the  upper  frequency  cutoff  of  the  laser  pulse  itself.  Assuming  the  first 
of  the  two  explanations,  and  assuming  a  vertical  exponential  distribution  of  the  source 
density  function,  the  1/e  depth  of  the  optoacoustic  source  is  estimated  to  be  0.084  mm. 
However,  since  the  length  of  the  laser  pulse,  at  3.5  ps,  has  a  cutoff  frequency  of 
around  300  kHz,  it  is  more  likely  that  this  is  the  cause  of  the  upper  frequency  cutoff. 
Currently,  we  have  no  way  of  determining  the  thickness  of  the  sound  source  layer  but 
the  above  results  show  that  it  must  be  equal  to  or  less  than  0.084  mm. 


Since  the  signal  was  bandpass  filtered  to  retain  only  the  frequency  components 
within  the  calibrated  frequency  range  of  the  hydrophones,  it  was  necessary  to  band 
pass  filter  the  theoretical  signal  over  the  same  band  before  direct  comparisons  could 
be  made.  The  E8  hydrophones  were  calibrated  from  100  kHz  to  1  MHz,  and  in  this 
band  their  response  is  known  to  be  flat.  Their  outputs  were  bandpass  filtered  from 
75  kHz  to  500  kHz  to  avoid  aliasing  because  the  data  recording  system  could  only 
sample  at  2  MHz.  The  H23  hydrophone  was  calibrated  from  1  kHz  to  150  kHz  and 
within  this  band  its  response  was  approximately  flat;  its  output  was  bandpass  filtered 
from  1  kHz  to  150  kHz. 

In  Fig.  14,  an  example  of  the  measured  acoustic  waveform  is  compared  with  the 
theoretical  model  prediction.  The  theoretical  model  includes  a  rectangular  bandpass 
filter  from  75  kHz  to  500  kHz  to  approximately  match  the  bandpass  filtering  of  the 
experimental  data.  The  prediction  was  based  on  a  laser  spot  size  of  8  mm2,  as 
measured  by  the  burn  method.  The  spot  size  inferred  from  the  high  cutoff  of  the 
spectrum  was  not  used  because  it  was  only  an  indirect  estimate  based  on  certain 
assumptions.  As  indicated  in  Section  II. D,  a  certain  amount  of  laser  energy  overhead 
loss  must  be  allowed  for.  There  was  no  way  to  measure  it,  but  it  was  found  by  trial  and 
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FIGURE  14  (a) 


COMPARISON  OF  AN  EXPERIMENTALLY  MEASURED  ACOUSTIC 
WAVEFORM  WITH  THE  CORRESPONDING  THEORETICAL  MODEL 
PREDICTION.  The  experimental  conditions  were: 
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error  that  a  value  of  0.13  J  gave  good  agreement  with  experimental  results.  This  figure 
is  expected  to  be  strongly  dependent  on  the  laser  spot  size.  In  Fig.  14(a),  the 
waveforms  are  compared.  There  is  much  ringing  and  this  is  mainly  due  to  the 
bandpass  filter.  The  theoretical  and  experimental  waveforms  appear  similar. 
However,  since  the  expected  acoustic  signal  is  a  short  impulse,  the  similarity  between 
the  theoretical  and  experimental  waveforms  is  mostly  a  product  of  the  bandpass 
filtering  rather  than  any  fundamental  agreement  between  theory  and  practice. 
Nevertheless,  the  comparison  is  encouraging  since  it  shows  that  the  model  predicts  a 
signal  pressure  level  that  is  in  approximate  agreement  with  the  experimental  result.  In 
Figs.  14(b)  and  (c),  the  corresponding  energy  density  spectra  are  compared.  The 
results  show  differences  in  the  shapes  of  the  spectra.  These  differences  are  most  likely 
due  to  the  oversimplicity  of  the  model. 
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FIGURE  14  (b) 

A  COMPARISON  OF  THE  ENERGY  DENSITY  SPECTRA 
OF  THE  THEORETICAL  AND  MEASURED  WAVEFORMS  IN  THE 
75  •  500  kHz  BAND  USING  THE  E8-59  HYDROPHONE 
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FIGURE  14  (c) 

A  COMPARISON  OF  THE  ENERGY  DENSITY  SPECTRA 
OF  THE  THEORETICAL  AND  MEASURED  WAVEFORMS  IN  THE 
1  150  kHz  BAND  USING  THE  H23  HYDROPHONE 
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IV.  CONCLUSIONS 


A  nonlinear  process  of  optoacoustic  sound  generation  was  modeled.  The 
process  was  one  of  inducing  an  explosive  reaction  at  the  water  surface.  The  explosive 
reaction  is  induced  by  using  a  high  intensity  laser  pulse.  The  experimental  results 
suggest  that  the  sequence  of  events  might  be  as  follows:  the  leading  edge  of  the  laser 
pulse  causes  a  phase  change,  most  probably  a  dielectric  breakdown,  within  a  thin 
surface  layer.  For  a  1.06  pm  laser  wavelength,  it  was  found  that  16  x  103  J/m2  was 
required  to  produce  the  phase  change.  As  a  result  of  the  breakdown,  the  remaining 
laser  energy  is  absorbed  by  the  layer.  Since  the  layer  is  very  thin,  the  absorbed 
energy  causes  an  explosive  expansion. 

A  blast  model  for  estimating  the  momentum  transfer  from  a  laser  pulse  to  a  solid 
surface12’13  was  adapted  to  give  an  approximate  model  of  the  optoacoustic  sound 
generation  process  in  water.  Bearing  in  mind  that  this  is  only  a  first  attempt,  the  model 
gave  estimates  of  the  acoustic  output  signal  which  were  in  reasonable  agreement  with 
our  experimental  results  as  well  as  some  of  those  of  Maccabee.1 1  The  signal  levels 
predicted  were  approximately  consistent  with  the  measured  data. 

There  were  a  number  of  discrepancies  between  theory  and  experiment.  The 
shape  of  the  signal  spectrum  as  predicted  by  the  theory  did  not  match  the  measured 
spectrum.  The  most  probable  cause  is  the  assumption  in  the  model  that  the  laser 
pulse  is  rectangular,  while  the  actual  laser  pulse  in  the  experiment  was  more  rounded 
with  a  short  rise  time  and  a  much  longer  fall  time.  The  measured  upper  cutoff 
frequency  in  the  off  axis  case  was  significantly  higher  than  the  corresponding  model 
prediction.  This  is  speculated  to  be  due  to  the  deformation  of  the  water  surface  which 
is  neglected  in  the  current  model. 

In  conclusion,  a  successful  first  attempt  has  been  made  to  model  a  nonlinear 
optoacoustic  sound  generation  process.  This  model  should  be  further  developed  to 
give  a  more  accurate  representation  of  the  process  and  thus  provide  a  better 
understanding  of  its  physics. 


V  FUTURE  PLANS 


Improvements  are  planned  in  both  the  theoretical  model  and  in  the  experimental 
apparatus 

The  theoretical  model  needs  to  be  improved  in  several  respects  The 
expenmental  results  indicated  that  a  portion  of  the  laser  energy  may  have  been  used 
to  produce  a  phase  change  within  a  thin  layer  of  the  water  before  a  blast  can  be 
generated;  this  needs  to  be  verified  by  a  quantitative  theoretical  analysis  The 
experimental  results  obtained  so  far  indicate  that  the  shape  of  the  laser  pulse  and  the 
deformation  of  the  water  surface  are  important;  therefore  the  model  should  be  modified 
to  take  these  into  account  Consequently,  a  model  for  the  deformation  of  the  water 
surface  is  also  required.  Photographic  data®  suggest  that  there  is  a  spherical 
indentation  of  the  surface  due  to  the  formation  of  a  gas  bubble  as  shown  in  Fig  15. 
Finally,  to  extend  the  validity  of  the  model  it  is  necessary  to  include  nonlinear 
propagation  effects  using  a  numencal  method  such  as  that  of  Hall  and  Holt.1®  Their 
method  uses  finite  element  analysis  in  a  dynamic  frame  of  concentric  elements,  as 
illustrated  in  Fig  16 

From  our  experience  in  the  laboratory,  a  number  of  improvements  in  the 
experimental  apparatus  are  required  in  order  to  provide  a  more  complete  set  of  data. 
An  optical  fiber  system  to  detect  laser  light  above  and  below  the  water  line  is  required 
to  monitor  *he  history  of  the  blast  and  verify  the  model  of  the  blast  process,  as 
illustrated  in  Fig  1  7  This  will  allow  us  to  determine  when  the  surface  has  reached  a 
state  of  dielectric  breakdown  It  may  also  allow  the  scattering  of  the  laser  light  by  the 
blast  to  be  studied  A  better  hydrophone,  with  a  wider  bandwidth  from  about  1  kHz  to 
10  MHz,  is  required  For  this  application,  commercially  available  shock  probes  may 
be  adequate  Finally,  to  keep  up  with  the  increased  data  rate,  the  data  acquisition 
hardware  needs  to  be  upgraded  to  handle  sampling  rates  of  at  least  10  MHz  per 
channel 

The  above  improvements  in  the  theoretical  model  and  the  experimental 
apparatus  will  enable  us  to  investigate  the  nonlinear  optoacoustic  processes  in  greater 
depth  We  will  be  able  to  observe  the  formation  of  the  blast  with  acoustic  and  optical 
sensors.  With  the  optical  sensors,  we  will  determine  the  proper  sequence  of  events  of 
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ANALYSIS  OF  NONLINEAR  PROPAGATION  IN  THE 
SPHERICALLY  SYMMETRIC  CASE 
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the  blast  process  and,  with  the  acoustic  sensors,  their  contributions  to  the  acoustic 
output.  The  process  is  expected  to  be  complicated  by  oscillations  of  the  gas  bubble, 
which  is  formed  as  a  by-product  of  the  blast.  When  the  process  is  understood  with 
sufficient  precision,  it  would  be  possible  to  make  estimates  of  its  optoacoustic  energy 
conversion  efficiency,  particularly  in  the  band  of  useful  sonar  frequencies. 


Calibration  charts  for  the  three  hydrophones  used  in  our  experiment,  the  E8-59, 
E8-62,  and  H23,  are  shown  in  Figs.  A-1  through  A-3.  The  H23  was  calibrated  for  the 
band  1  kHz  to  150  kHz,  the  E8-59  was  calibrated  from  100  kHz  to  1  MHz,  and  the 
E8-62  from  200  kHz  to  1  MHz.  These  charts  were  made  by  Naval  Research 
Laboratory,  Orlando,  Floriua.  Although  the  calibration  chart  for  the  E8-62  only  extends 
to  200  kHz,  we  were  able  to  estimate  its  calibration  down  to  100  kHz  by  comparing  the 
broadband  response  of  the  E8-62  to  that  of  the  E8-59.  Between  200  kHz  and  600  kHz 
the  responses  of  these  two  hydrophones  are  almost  the  same,  and  therefore  we  were 
able  to  get  a  relative  calibration  of  the  E8-62  by  taking  the  Fourier  transforms  of  their 
responses  to  two  identical  broadband  signals  which  had  the  same  levels  beyond 
200  kHz  and  then  dividing  the  response  of  the  E8-62  response  by  that  of  the  E8-59  for 
frequencies  between  100  kHz  and  200  kHz.  The  hydrophones  are  shown  in  Figs.  A-4 
and  A-5.  The  dotted  regions  represent  the  active  elements. 
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FIGURE  A-1 

CALIBRATION  CURVE  OF  THE  E8-59  HYDROPHONE 
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FIGURE  A-2 

CALIBRATION  CURVE  OF  THE  E8-62  HYDROPHONE 
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FIGURE  A-4 


THE  E8  HYDROPHONE:  DIMENSIONS  IN  CENTIMETERS 
The  dotted  area  indicates  the  acoustically  active  region 
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The  Q-switch  driver  control  circuit  is  a  TTL  circuit  designed  to  send  a  series  of 
trigger  pulses  to  the  Q-switch  driver.  It  can  be  set  to  send  a  designated  number  of 
pulses  with  a  particular  repetition  rate  after  a  preset  delay.  The  control  circuit  has  four 
sections  to  it  as  shown  in  Fig.  B-1.  The  first  section  buffers  the  trigger  pulse  from  the 
laser  firing  unit  and  provides  a  clean  trigger  pulse  to  the  state  graph  controller.  This 
section  is  necessary  because  the  laser  will  often  send  two  or  more  trigger  pulses.  The 
second  section  is  the  state  graph  which  controls  the  sequencing  of  events.  The  third 
section  is  the  timing  section,  which  defines  the  number  of  pulses,  the  repetition  rate, 
and  the  initial  delay.  The  fourth  section,  a  clock,  provides  the  timing  for  the  state  graph. 
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FIGURE  B-  I 

SCHEMATIC  OF  THE  Q-SWITCH  DRIVER  CONTROL  CIRCUIT 
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APPENDIX  C 

DETERMINATION  OF  SOURCE  WIDTH  AND  DEPTH 


A!  any  given  frequency,  the  discrepancy  between  the  acoustic  signal  level 
omc.ced  by  a  distributed  source  and  a  point  source  of  equal  power  is  termed  a 
q  •"  k*  cm  css  “  For  a  white  noise  source,  diffraction  losses,  resulting  from  the  finite 
i  "  >-■  s  rs  o*  t^e  source  volume,  will  impose  an  upper  cutoff  frequency  on  the 
’  c  s  gna:  At  frequencies  above  the  cutoff,  significant  diffraction  loss  is  incurred 
t  •  ,He  destructive  interference  of  signal  elements  from  spatially  separate  points 
a . e  sc..'ce  vO'ume  Conversely  if  the  cutoff  frequency  is  observable  then  it  is 

•  -  ^a«e  es*  ^afes  of  the  physical  dimensions  o<  the  source  volume  The 

•  "  ,-s  a’ *~e  *'equency  w^ere  the  acoustic  wavelength  is  approximately  equal  to 

"  ••  .  *-  *ed  e-gh'  o4  *he  source  >n  the  direction  of  the  observer  The  exact  cutoff 

.■r  f  e  soecva.  response  s  dependent  on  the  details  of  the  source  density 
a  ••  ”'e  snu'ce  vo  sme  Nevertheless  it  is  possible  to  get  an  approximate 

•  •  ■  ••  *  sc., me  /O'ume  dimensions  by  approximating  the  source  density  with  an 

•  -*  c'da-d  ‘.mr’cn  SuC"  as  a  Gaussian  or  an  exponential  function  The 
■  i  -u'Awen  sou'ce  d  versions  and  the  cutoff  frequency  for  a  few 

■  u  i«’'s  4  ,t'  rns  are  derived  below 


source 


A.  i  he  uniform  Density  f-unction 

Assuming  a  uniform  surface  source  density  function,  the  diffraction  loss  factor  Da 
is  given  by 


D- 


' 


/*  x  =  r  _ 

(nr2)'1  2V(r2-x2)  cos  ((x/X)  27rsin0Q)  dx 

x  =  -r 


(C.1) 


where  r  is  the  radius  of  the  blast  aperture,  80  is  the  angle  to  the  receiver  referenced  to 
the  downward  vertical,  and  X  is  the  wavelength  corresponding  to  the  cutoff  frequency. 
Using  the  substitution,  u  =  x/r,  du  =  d xJr,  Eq.  (C.1)  becomes 


Da  =  J  (nr2)"1  2rV(  1  -  u2)  cos  ((ur/X)  2resin0o  ) 


(rdu) 


(C.2) 


Simplifying  Eq.  (C.2)  and  taking  advantage  of  the  function’s  symmetry  about  the  origin 
gives 


Da 


-1  _ 

(4/k  )V(  1  -  u2)  cos  (u2nrsin0o/X)  du 

J0 


The  solution  is  expressible  as  a  Bessel  function, 


(C.3) 


Da  =  (  X/(nrsin0o) )  J-|(27trsin0o/X)  ,  (C.4) 

where  is  a  first  order  Bessel  function.  By  substituting  the  wavelength  X  with  the 
-3dB  cutoff  wavelength  ^  which  is  related  to  the  cutoff  frequency  fc  by 

Xc  =  c/(2rcfc)  ,  (C  5) 

where  c  is  the  speed  of  sound  in  water,  and  setting  Da  equal  to  (1/V2),  an  expression 
for  the  radius  r  of  may  be  found. 


59 


A  Gaussian  function  is  another  plausible  approximation  for  the  surface  source 
distribution.  The  equation  is 


Da  -(wi*) 


2\-1 


o 

exp[*(x2  +  y2)/a2]  cos(x2x  sin0Q  Ik)  dx  dy 


(C.6) 


where  a  is  the  1/e  radius  of  the  blast  aperture.  The  integration  with  respect  to  x  gives 
Oa  -  2  ( W  /"  «pHa,  Sine,  oft  I  dy  <0.7, 

Integrating  with  respect  to  y  and  simplifying  gives 

Da  =  exp[-(an  sin0o  /X.)2]  (C.8) 

Setting  Da  equal  to  1/V2,  X  to  kc  and  solving  for  a  gives 

a  *  V(ln2/2)  ^/(nsm©^  (C.9) 


C  The  Exponential  Density  Function 

The  exponential  density  model  is  a  plausible  approximation  for  the  vertical 
source  density.  The  diffraction  loss  for  the  0o=O  case  is  given  by 

r * 


exp[-z/a]  exp(j2nz/X]  dz 


(0.10) 


'  0 


The  solution  is  simply 


Setting  Da  equal  to  1W2,  Xto  5^,  and  solving  for  a  gives 


a-X/2*  (C.l  2) 

D.  Estimates  of  Width  and  Depth  of  the  Ootoacoustic  Source 

Estimates  of  width  and  depth  estimates  based  on  data  at  90  =  45°  and  90  =  0° 
obtained  from  the  three  different  models  of  source  density  are  given  in  Table  C-1 . 


TABLE  C-1 

ESTIMATED  DIMENSIONS  OF  OPTOACOUSTIC  SOURCE 


density  function 


width 
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